Turbulence lifetimes: what can we learn from the physics of glasses? 
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In this note, we critically discuss the issue of the possible finiteness of the turbulence lifetime 
in subcritical transition to turbulence in shear flows, which attracted a lot of interest recently. 
We briefly review recent experimental and numerical results, as well as theoretical proposals, and 
compare the difficulties arising in assessing this issue in subcritical shear flow with that encountered 
in the study of the glass transition. In order to go beyond the purely methodological similarities, we 
further elaborate on this analogy and propose a qualitative mapping between these two apparently 
unrelated situations, which could possibly foster new directions of research in subcritical shear flows. 

PACS numbers: 47.27.Cn, 47.27.eb, 64.70.P- 



I. INTRODUCTION 

Statistical physics has devoted a lot of effort to the 
study of fully developed turbulence, but much less to that 
of the transition to turbulence pQ , which occurs when the 
Reynolds number, the ratio of the advection time to the 
viscous time, is increased. The transition is commonly 
observed in flow regimes lacking linear instability and is 
referred to as globally subcritical [2HI]- 

The plane Couette flow, driven by two plates moving 
parallel to each other in opposite directions is linearly 
stable at all Reynolds number and, as such, is the epit- 
ome of globally subcritical transitions [5]. Other flows 
usually transit to turbulence before linear instability sets 
in. These include the circular Poiseuille flow (cPf) and 
the plane Poiseuille flow (pPf), which are driven by a 
pressure gradient respectively along a circular pipe or be- 
tween two parallel plates, as well as the counter-rotating 
Taylor Couette flow (TCf), driven by two concentric 
cylinders rotating in opposite directions. In all these 
cases, the transition is particularly delicate to understand 
owing to its abrupt character. Many non-trivial unsta- 
ble solutions of the Navier-Stokes equations coexist in 
phase space with the laminar flow, leading to a complex 
spatio-temporal dynamics in real space including the nu- 
cleation and the growth or decay of turbulent domains 
called 'puffs' (pPf) or 'spots' (pCf) -see, e.g., U\ for 
cPf, [ErEa for pCf, [IT] for pPf, and [HE2 f »r T Cf- 

A recent surge of interest has been motivated by the 
audacious proposal that shear flow turbulence could re- 
main transient up to arbitrarily large Reynolds number, 
opening ways towards a better control of such turbu- 
lent regimes |14j . This proposal was motivated by new 
experimental and numerical observations in cPf |14j re- 
garding the statistics of turbulent lifetimes, in contradic- 
tion with those previously obtained in cPf [fSJ [TB] and 
pCf jTTl [18] . This contradiction, together with the exper- 
imental discovery of a spectacular long wavelength peri- 
odic organization of the laminar-turbulent coexistence in 
pCf and TCf [19l [20], has motivated new experiments 
in TCf [21] and cPf [221 [23], the development of various 
models [24-28 and an impressive number of numerical 



studies [2"2l f^TTfTTU] . As a result, a better knowledge of the 
organization of the phase space and some comprehension 
of the mechanisms at play in the coexistence dynamics 
has been gained. However, whether the turbulence life- 
time remains finite or diverges at some critical Reynolds 
number is still a matter of debate. 

Interestingly, the possible existence of finite, yet ex- 
tremely large relaxation times, together with the coexis- 
tence in phase space of many unstable solutions (saddles 
from a geometric point of view), is reminiscent of the 
physics of glasses (see, e.g., [1TH4*3] ). Whether the struc- 
tural relaxation times of a glass really diverges at a given 
finite temperature, or remains very large but finite at any 
positive temperature, is an important issue related to the 
existence of a genuine phase transition to an ideal glass 
state. After several decades of intense activity, a major- 
ity of the community has come to the general agreement 
that no fitting procedure, however cautious, could bring 
a definitive answer to that question [U] . In spite of this 
apparent failure, new (and perhaps even more interest- 
ing) questions have appeared as a result of this research 
effort, and have driven the field of glasses towards impor- 
tant conceptual progresses |4"5] . 

Here, in the legacy of Y. Pomeau [JS], we would like to 
draw again the attention of the statistical physics com- 
munity towards the old standing problem of the transi- 
tion to turbulence. On the one hand, many new concepts 
have been developed in the physics of glasses (and more 
generally of disordered systems) in the last twenty-five 
years, since the publication of this seminal paper. On the 
other hand, the development of experimental techniques 
such as Particle Image Velocimetry, and the exponential 
increase of the numerical capacities now give access to 
never achieved before datasets about the transition to 
turbulence. In such circumstances, one may hope signif- 
icant progress to emerge again from the confrontation of 
both fields. 

As a first modest contribution in this direction, we wish 
to discuss, in the light of a simple model inspired by the 
so-called Random Energy Model [46[ |47] introduced to 
describe disordered systems, what scenarios could lead 
to diverging turbulence lifetimes at some finite Reynolds 
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number. Our aim here is to illustrate the possible rich- 
ness of the analogy, and to possibly inspire new exper- 
imental protocols and new data analysis methods. It 
might also point at deeper connections between these two 
active fields of nonlinear and statistical physics. The pa- 
per is organized as follows. We shall first propose a con- 
cise review of the major results dealing with subcritical 
transition to turbulence. We then introduce our model, 
before discussing the possible outcomes of the analogy. 

II. FINITE VERSUS DIVERGING LIFETIMES 

We first briefly review the experiments and direct nu- 
merical simulations reporting the increase of the turbu- 
lence lifetime when the Reynolds number is increased. 
We then describe two complementary numerical ap- 
proaches, a low dimensional (mean field-like) one, focus- 
ing on the phase space structure, and a spatially extended 
one focusing on the real space dynamics. 

A. Experimental and Numerical observations 

To our knowledge, the first systematic measurement of 
turbulence lifetimes was conducted in the pCf [T71 [TS] . 
Two different kinds of experiments have been performed, 
differing by the way the initial condition is prepared. 
In what we shall call type-A experiments, the Reynolds 
number, i?, is set to its value of interest and the laminar 
flow is disturbed locally at the initial time. In type-B 
experiments, a turbulent flow is prepared at high i?, and 
quenched at the initial time down to the R value of inter- 
est. In both cases, one monitors the evolution of the tur- 
bulent fraction frit), which characterizes the coexistence 
dynamics of laminar and turbulent domains (see Fig [lj . 
For R> R g , fr(t) fluctuates around some average value, 
which remains finite on experimental timescales. For 
R < R u , fr{t) relaxes towards zero, without display- 
ing any transient regime. In between, for R u < R < R g , 
frit) exhibits a first rapid decay, followed by a long tran- 
sient quasi-steady regime, before a large fluctuation sets 
it to zero. The lifetime of these transients are exponen- 
tially distributed and the average value r was reported 
to diverge like (R g — iJ) -1 . 

The cPf was later investigated in various ways. In [6j 
116] a puff was generated inside a constant flow rate pipe 
flow by introducing a short duration perturbation. Then 
R was reduced and the subsequent evolution of the puff 
was monitored as it progressed downstream. The proba- 
bility of observing a localized disturbed region of flow as 
a function of distance downstream is exponential and the 
time required for half of the initial states to decay, ti/ 2 , 
was reported to diverge like {R c — R)~ x , in agreement 
with the observations made in the pCf. Other protocols 
lead to the same conclusions [HJ [TC] . 

However, these results were challenged later by an- 
other experimental study [14] . In a pressure driven flow 




FIG. 1: Typical snapshots of the laminar-turbulence coexis- 
tence at intermediate Reynolds number in plane Couette flow. 
Top: turbulent spot in a small aspect ratio setup, created by 
a localized perturbation. Bottom: Large scale coexistence 
in a large aspect ratio setup, following a quench from high 
Reynolds number. 

through a very long pipe, the authors could record much 
longer dimensionless observation times. They could de- 
termine the probability to be turbulent after a time pe- 
riod given by the distance between the perturbation loca- 
tion and the outlet, as a function of flow rate. For short 
times, the data are within the error bars of [fH [16] but 
for longer times they deviate from the divergent behavior 
reported above and are better represented by an expo- 
nential variation: r = cxp(ai? + 6), without singularity 
(here and in what follows, a and b denote generic fit pa- 
rameters) . Finally in a recent experimental study of tur- 
bulence in pipe flow spanning height orders of magnitude 
in time, drastically extending all previous investigations, 
it was claimed that the turbulent state remains transient, 
with a mean lifetime, which depends super-exponentially 
with the Reynolds number: r oc exp(exp(ai? + b)) [48 . 

Intense numerical simulations of the cPf have also been 
conducted, but did not clarify the situation. In [HOES] a 
diverging behavior of the turbulent lifetimes compatible 
with the experimental results of [6, 16 is reported. Later 
in [31] . the authors -one of which is common to |15|- 
conducted new simulations and reanalyzed older data, 
concluding to an exponential dependence such as the one 
reported in 14 . Altogether despite intense experimental 
and numerical effort, no definitive answer regarding the 
divergence or finiteness of turbulence lifetime could be 
obtained from the fit of data by phenomenological func- 



tional forms. 

As stated in the introduction, this issue is not spe- 
cific to the transition to turbulence. When a liquid is 
suddenly quenched below its crystallization temperature 
and if crystallization can be avoided, the liquid enters a 
state, called supercooled liquid, in which the relaxation 
time increases by several orders of magnitude over a lim- 
ited range of temperature [H]. A divergence at a finite 
temperature of the relaxation time would signal an ideal 
glass transition, and would thus be of high interest, at 
least at a conceptual level. Despite huge efforts made 
to measure the variations of the relaxation time over an 
experimental window of more than ten decades, no clear 
consensus has been obtained yet. More precisely, the 
available data are both consistent with fits including a 
divergence at a finite temperature T c > 0, and with fits 
diverging only at T = 44 . 

The same difficulty is also expected to appear in the 
context of turbulence. We illustrate this point on experi- 
mental data recently obtained in the case of the TCf [3T] , 
when only the external cylinder is rotating. The TCf is 
then, like the pCf, linearly stable for all R. Also, because 
the TCf is a closed flow, one can record very long times. 
In this experiment, the angular rotation of the external 
cylinders fixes the Reynolds number. The flow was per- 
turbed by rapidly accelerating the inner cylinder in the 
direction opposite to the rotation of the outer cylinder 
and immediately stopping it. After a short regime of 
featureless turbulence, the flow exhibits long transients 
characterized by the coexistence of laminar and turbulent 
domains, before eventually relaxing towards the laminar 
flow. The distribution of lifetimes is again exponential, 
and the authors argue that the mean turbulent lifetime 
does not diverge and rather behaves in the transitional 
regime as a double exponential r oc exp(exp(ai? + b)), as 
observed in the cPf [4"5] , 

It is interesting to note that in the oldest experiments, 
the debate about the functional dependence of the aver- 
age turbulent lifetime on the Reynolds number was con- 
centrating on the choice between the two following forms: 



r/r = exp(R/Ro) 
, R 

r/T 



Rr ~ R 



a > 0, 



(1) 
(2) 



whereas the most recent experiments, both in the case 
of the cPf [J5] and the TCf [3T], have access to much 
longer experimental timescales and point at a double ex- 
ponential behavior. This last functional form ensures 
a very fast increase of r without singularity, and could 
give the impression that it solves the above debate. How- 
ever as learnt from the physics of glasses, the debate has 
actually been shifted towards two alternative functional 
forms, namely: 



ln(r/r ) = \exp(R/R 



(3) 

a > 0. (4) 
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FIG. 2: Probing finite lifetime experimentally (color online): 
Relaxation lifetimes of turbulent initial conditions in a Taylor- 
Couette flow, with rotating external cylinder and internal 
cylinder at rest (data from [21]). Four possible fits are pro- 
posed as indicated in legend. Top and bottom panels are in 
lin-lin and lin-log scales respectively. Times are given in units 
of d/roujo, where d is the gap between the two cylinders, ro is 
the radius of the external cylinder and u>o its angular veloc- 
ity. All fits were performed using a standard least square fit 
procedure. Fit with Eq. (1) was obtained imposing a linear 
fit of log(r) vs. R. Fit with Eq. (2) was obtained imposing 
a linear fit of log(r) vs. log(R c /(R c — R)), optimizing the fit 
quality on R c . Fit with Eq. (3) was obtained imposing a lin- 
ear fit on log(r) vs. exp(R/Ro), optimizing the fit quality on 
i?o. Finally, fit with Eq. (4) was obtained imposing a linear 
fit of log(r) vs. R C /(R C — R), optimizing the fit quality on 
R c . In this last case, a was thus set to one. It was checked 
that other values of a up to 3 cannot be discriminated. For 
Eq. (1), the regression coefficient r : 2 is equal to 0.971. For all 
the other cases, ri = 0.994. 

As a matter of fact, Eq. Q, which has (to the best of our 
knowledge) not yet been proposed in the context of the 
transition to turbulence, is a very standard form called 
the Vogel-Fulcher-Tammann (VFT) law in the physics of 
glasses [33- 
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Figure [2] displays the data obtained in TCf [UJ -which 
are available online as supplementary material- together 
with possible fits by the four functional forms proposed 
above. Note that we have only performed global fits of 
the data, without trying to extract various regimes and 
crossovers as can be done in the case of glasses [44] . One 
clearly observes that indeed the relevant variable to de- 
scribe the growth of the turbulent lifetimes is hi(r), as 
soon as really large times are considered. However, one 
also sees that apart from the simple exponential form 
Eq. [T] all other descriptions are not discriminablc, so 
that there is no definitive way to rule out or confirm the 
existence of a singularity. One faces the same difficulty 
as in the physics of glasses: the lifetimes to be measured 
become very large, which makes it difficult to accumu- 
late significant statistics. The experimental results are 
thereby confined in a finite range of Reynolds number 
or temperature, from which even with high-quality data 
sets, the divergence of a characteristic time cannot con- 
vincingly be determined from fits. 

Before concluding this section, let us mention that the 
double-exponential form Eq. [3] has been justified on the 
basis of extreme value statistics [49 . However, as stated 
by the authors, the argument is only local, as it involves 
an expansion in R around a given reference value. Hence 
no clear conclusion can be drawn from the theoretical 
argument of [49] on the issue of the divergence of r at 
a finite or infinite value of R. Finally, let us empha- 
size that for now, we have left aside all issues related to 
finite size effects, which in turn can severely alter the 
functional dependence of time and length scales in tran- 
sitional regimes. 



B. Low dimensional vs. spatially extended models 

The globally subcritical transition to turbulence is by 
definition controlled by solutions of the Navier-Stokes 
equation, which do not branch continuously from the 
laminar flow solution when the Reynolds number is in- 
creased. These solutions -of various kinds, station- 
ary states, traveling waves, or more complex coherent 
structures- are unstable and form hyperbolic states, with 
stable and unstable manifolds. Early indication of the ex- 
istence of these solutions were reported in pCf, both nu- 
merically 50 52] and experimentally [53] , The intricate 
network made of these manifolds and their connections 
then serves as a skeleton for the turbulent flow. 

In principle one would like to collect all such states, 
estimate their dynamical weight and calculate statistical 
averages from periodic orbit theory |54j . In practice, one 
must restrict the analysis to low-dimensional models (551 — 
159] or to simulations [H HS1 [Ml El ES [HO] , performed 
in the so-called minimal flow unit assumption [ST]. Do- 
ing so, it was shown that the regions of initial conditions 
for which long lifetimes exhibit strong fluctuations and a 
sensitive dependence on initial conditions were separated 
from the regions with short lifetimes and smooth varia- 



tions by a border, the so-called "edge of chaos". Later, 
some exact solutions with codimension-1 stable mani- 
folds have been identified as edge states, that is, solutions 
that locally form the stability boundary between laminar 
and turbulent dynamics |34H36| . These important results 
contributed to make concrete the picture borrowed from 
dynamical system theory of a turbulent repellor, sepa- 
rated from the laminar state by a set of edge states con- 
nected through heteroclinic manifolds. In particular the 
existence of the above non trivial solutions has served 
to understand the exponential distribution of lifetimes in 
the transitional regime. 

Unfortunately the above picture does not bring a 
definitive answer to the issue of the existence of a fi- 
nite critical Reynolds number above which turbulence is 
sustained. As argued in J35J 123 E2] , the reason is that 
the dynamics being either projected on a small set of 
modes or limited to small computational domains with 
periodic boundary conditions, it cannot capture the gen- 
uinely spatiotemporal coexistence of laminar and turbu- 
lent states observed in open and unbounded flows. In 
particular, it can neither capture the long wavelength 
modulation of turbulent intensity, nor the regime of al- 
ternating laminar and turbulent stripes, first observed 
experimentally in pCf and TCf [JJ|1 [3U] then reproduced 
numerically in pCf [29j E2 ESI EH1 HOj [63] . 

As a matter of fact it is known for long that, accord- 
ing to the scenario called spatiotemporal intermittency 
(STI [HI]), transient chaotic local states of a distributed 
system may evolve into a sustained turbulent global state 
due to spatial couplings [65H68] , Following this path, it 
was demonstrated that a simple ID-model of cPf, com- 
posed of coupled maps, does indeed captures remark- 
ably well the character of the turbulent pipe flow in the 
transitional regime and exhibits a critical transition to- 
wards sustained turbulence via spatiotemporal intermit- 
tency [28 . The transition is further believed to belong 
to the Directed Percolation class [651 E3 [69] , as already 
suggested in [15] for pCf, and recently reconsidered in 
cPf [70]. 

Finally, it was shown by means of fully resolved direct 
numerical simulations of the Navier Stokes equation, that 
there exists a crossover length scale of the order of 10 2 
times the cross-stream length below which the spatio- 
temporal processes at play in large-scale simulations and 
experiments are not captured |27] . Since then, a number 
of numerical investigations of large aspect ratios cPf and 
pCf have reproduced the complex spatio-temporal coex- 
istence of laminar and turbulent states, and identified the 
first hydrodynamics mechanisms at play [22 EH EH SO] • 



C. Summary 

Let us now recast the essential features of subcritical 
transition to turbulence, which will feed the analogy with 
the glass transition that we would like to propose and 
illustrate on a simple model: 



5 



• Subcriticality: While the laminar flow is stable 
against infinitesimal perturbations, finite ampli- 
tude perturbations may trigger an abrupt transi- 
tion towards a disordered flow. Such disordered 
flow can also be obtained by quenching fully tur- 
bulent flows. 

• Spatio-temporal intermittency: This disordered 
flow is made of turbulent domains, which move, 
grow, decay, split and merge leading to spatio- 
temporal intermittency, that is a coexistence dy- 
namics in which active/turbulent regions may in- 
vade absorbing/laminar ones, where turbulence 
cannot emerge spontaneously. 

• Transients and M eta- stability: For large enough 
Reynolds number this disordered flow has long 
lifetimes, which are distributed exponentially. 
Whether the associated characteristic time diverges 
at a finite Reynolds number is still a matter of de- 
bate. For low Reynolds number, say R < R u , or 
small enough disturbances, the flow relaxes rapidly 
towards the laminar flow. 

• Unstable states: When increasing the Reynolds 
number a larger and larger number of unstable 
finite amplitude solutions appear in phase space. 
Some have been identified as edge states separat- 
ing the others from the laminar state. 

On the basis of these observations, we now propose an 
analogy with the glass transition, the goal of which is 
twofold. It offers a new framework to discuss possible 
scenarios which could lead to diverging turbulence life- 
times at some finite Reynolds. In addition, it may also 
inspire new analysis of the subcritical transition to tur- 
bulence, and might drive this field of research towards 
important, if not conceptual, methodological progresses. 

III. A SIMPLE MODEL FOR THE 
DIVERGENCE OF LIFETIMES 

A first analogy with the study of glasses has been em- 
phasized above at the methodological level, when study- 
ing the divergence of turbulence lifetime from empirical 
data. In this section, we further develop the analogy, and 
show that it can also be fruitful to some extent at the con- 
ceptual level, within the complex landscape picture. For 
illustration purposes, we introduce a simple model that 
captures, as an essential ingredient, the wandering on a 
complex landscape, and yields interesting predictions for 
the dependence of the turbulent lifetime on the Reynolds 
number. 



A. The landscape picture 

One of the main interests of the analogy with glasses is 
related to the energy landscape picture [7T]. In glasses, 



the transition between the liquid state and the crys- 
tal occurs, when lowering the temperature, through the 
emergence of a slowly relaxing glassy state (the super- 
cooled liquid) [41] . resulting from the wandering of the 
phase-space point representing the system in a complex 
energy landscape (72], mostly composed of many sad- 
dles (321 [721 [73] (though local minima also play an impor- 
tant role at low enough temperature). The most striking 
feature of the glass transition, the rapid increase of the 
relaxation time, by several orders of magnitude over a 
moderate range of temperature is interpreted as a conse- 
quence of this complex dynamics in phase space. 

Given the experimental and numerical evidences that 
we have reviewed in the previous section, it is tempting to 
see in the complex phase-space of the spatially extended 
dynamical systems an analog of the complex energy land- 
scape of glasses. Though numerical investigations of tur- 
bulent flows have not been able yet to characterize the 
number of unstable solutions as a function of the volume 
of the flow, the analogy with glasses suggests that this 
number of solutions may grow exponentially with the vol- 
ume of the system. Characterizing the state of the flow 
by its turbulent energy per unit volume, e = E/V (that 
is the excess kinetic energy with respect to the laminar 
flow), we assume that the number f2y(e, R) of unstable 
solutions at a given energy density e and Reynolds num- 
ber R grows exponentially with the volume V according 
to 

tl v (e,R)~e v '^ R \ (5) 

thus defining an entropy density s(e, R). At low Reynolds 
number, no unstable states exist, so that we assume that 
the entropy s(e, R) is equal to zero for all e > if R is less 
than a characteristic value R u . For R > R u , we assume 
that s(e,R) > on an interval e m i n (R) < e < e max (R), 
and s(e,i?) = otherwise, meaning that unstable states 
exist only in the energy interval (£ m i„,£ mlB: ). 

B. Diffusion in the energy landscape: the REM 
analogy 

Turning to dynamics, we assume on the basis of the 
experimental and numerical observations reported in sec- 
tion |II B| that the turbulent flow spends most of its time 
close to unstable solutions, and that the evolution of the 
flow can be considered as a succession of jumps between 
different unstable solutions. If however the flow ends up 
in the laminar state, the evolution stops until an external 
perturbation is imposed. Taking into account the pres- 
ence of the absorbing laminar state is obviously essential 
to determine the lifetime of the turbulent flow. This will 
be the focus of Sect. |III C| Yet, in a first stage, it is in- 
teresting to consider the evolution of the turbulent flow 
in the absence of the absorbing laminar state, in order to 
make the analogy with glass models emerge more clearly. 

As it is unlikely that a large amount of energy could 
be injected or dissipated within a short time period, one 
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expects that the energy of successively visited unstable 
solutions are close to one another. At a coarse-grained 
level, it is then natural to assume that the energy e 
evolves diffusively. In order to take into account the 
variation with e of the number of unstable states, the 
evolution should also be biased toward values of the en- 
ergy with a high entropy s(e,R). More precisely, the 
bias should depend on the derivative of the entropy with 
respect to the energy (a constant entropy introduces no 
bias in the dynamics). Altogether, the simplest evolu- 
tion equation for the energy e(t) incorporating the above 
ingredients is the following Langevin equation 



ds 
~dt 



(6) 



where the prime denotes a derivative with respect to e. 
To enforce the finite range of values £ m i n < e < £ max , 
reflecting boundary conditions are assumed at e m i n and 
£max- The parameter 7 is a proportionality coefficient to 
be determined later on, included for dimensional reasons. 
The term A accounts for dissipative effects, and £(t) is a 
white noise describing the energy injection mechanism, 
satisfying 



(7) 



where D is a diffusion coefficient in energy space. These 
are obviously strong simplifications: the dissipation rate 
could in principle depend on e and the noise should rather 
be considered as colored and multiplicative in such non- 
equilibrium systems, but we wish to keep the model as 
simple as possible for the sake of illustration. The as- 
sumption of a constant dissipation rate is however jus- 
tified in the limit where the width £ max — £ m in of the 
accessible energy range is small with respect to £ m in- Be- 
sides, considering that the noise is self-generated by the 
turbulent fluctuations, and thus results from the super- 
position of a number of contributions proportional to the 
volume V, one expects the diffusion coefficient to scale as 
D = Dq/V. Note that all parameters 7, A and Dq may 
depend on the Reynolds number R. 

The Fokker-Planck equation describing the evolution 
of the probability distribution P{e, R, t) then reads 



dP 

~dt 



de ' 



;{hs'-\]P 
The stationary solution P(e, R) is obtained as 



D d 2 P 



P(s, R) = \ cxp (J^ ( 7S (e, R) - Xe) 



(8) 



(9) 



Following standard statistical physics arguments, one ex- 
pects the distribution P(e,R) to be proportional to the 
number of unstable states fV(e, R) ~ e^' 6 '^, which 
imposes 7 = D a . Introducing the parameter (3 = A/Do, 
the stationary distribution then reads 



1 
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FIG. 3: (Color online) Sketch of the entropy surface and its 
slopes along the energy density direction, together with the 
path followed by the flow while varying the reynolds num- 
ber. Colors of the surface go from blue to red with increasing 
Reynolds number. The green line indicates e m i n (R). Fixing 
some Reynol ds number R, a given dark line, one sets P(R). 
Solving Eq. (Ill then graphically amounts to finding a slope 
along the energy density direction equal to /3(R). A solu- 
tion exists if P(R) < Po(R), the slope at the intersect with 
Emin(R)- Varying R, one follows the blue path on the surface, 
eventually leading to the value R g such that fi(R g ) = /3o(R g ). 



where we have emphasized the independence of the pa- 
rameter fi, which in the present context describes the 
balance between the energy injection and the dissipa- 
tion, as does the inverse temperature at equilibrium. 
If V is large, the distribution is dominated by the en- 
ergy e(R) which maximizes the argument of the expo- 
nential, namely s(e, R) — j3(R)e. If the maximum of 
s(e, R) — f3(R)e lies within the interval £ m i n (-R) < £ < 
£max(-R), the most probable energy is the solution of 



s'{e(R),R) - f3{R) = 0. 



(11) 



P(e, R) = - exp[F( S (e, R) - (3(R)e)] 



(10) 



Assuming the entropy s(e, R) to be a concave function 
of e (see figure [3]), s'(s,R) is a decreasing function of e, 
and thus has its maximum at £ = s m - ln (R). 

We now introduce the key element of the model, which 
we borrow from the Random Energy Model (REM) [46, 
I47j . one of the simplest disordered models exhibiting a 
freezing transition. The specificity of the REM, which is 
the ingredient leading to a glass transition in this model, 
is that the entropy has a finite slope at the minimum 
energy. We shall now see what are the consequences of 
this freezing (or glass) transition in the present model for 
turbulent flows. 
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By analogy, we thus assume that s'(e,R) takes a fi- 
nite value, denoted as Po{R), when e — > e m i n (R). As 
a result, if /3(R) < /3 {R), Eq. (Ill generally admits a 
solution e(R) > £ m in(-R)- In contrast, if (3(R) > (3o{R), 
Eq. (11) has no solution, and s(e,R) — f3(R)s is maxi- 
mum at e = e m j n (i?). The probability distribution then 
concentrates on e m - m (R). Intuitively, one expects (3(R) 
to be a decreasing function of R (that is, the tempera- 
ture characterizing the fluctuations increase with the 
Reynolds number). On the other hand, the total number 
of unstable states increases with the Reynolds number, 
and it is thus plausible that (3o(R) increases (or at least 
remains constant) with R. This suggests the existence of 
a Reynolds number R g such that (3(R g ) = (3o(R g ). In this 
case, the average energy e(R) is larger than £ m j n (i?) for 
R > R g , while the dynamics in phase space concentrates 
on the states of minimal energy for R < R g . 

As emphasized at the beginning of this section, these 
conclusions hold under the unphysical hypothesis that 
no laminar state is present. However, if the paths lead- 
ing from the unstable states to the laminar one are rare 
enough, the flow is likely to visit a large number of unsta- 
ble states, and should thus partially equilibrate, before 
ending up into the laminar state. It is then plausible 



that the equilibrium distribution given in Eq. ( 10 ) qual- 



itatively describes this quasi-equilibrium regime. A nat- 
ural assumption is that most of the paths leading to the 
laminar state are connected to unstable states close to 
£min(R), the so-called edge states in the context of tur- 
bulence. As for R < R g , the average energy remains 
close to e m i n (R), the flow should reach the laminar state 
in a reasonably short time. Conversely, for R > R g , the 
typical energy remains well above the threshold e m - ln (R) , 
and one expects that it takes a very large time to find 
the laminar state, as it implies excursions very far from 
the typical energy. 

Hence, the Reynolds number R g appears as a transi- 
tion (or crossover) value between a regime of short tur- 
bulent lifetime and a regime of large lifetime. Note also 
that the turbulent lifetime should essentially vanish be- 
low the Reynolds value R u where unstable states cease 
to exist. 



C. Divergence of turbulent lifetime 

In this section, we now try to put the above argu- 
ments on a more quantitative basis. We define the tur- 
bulent lifetime as the mean time to reach the laminar 
state after a sudden quench from a higher Reynolds num- 
ber value, where turbulence is established. This situation 
can be modeled using Eq. ([6| for the stochastic dynam- 
ics of e(t), with now an absorbing (instead of reflecting) 
boundary at e = e m ; n to account for the presence of the 
laminar state. The initial condition at t = is chosen 
as e(0) = £ maX j to model the quench from high energy 
turbulent states. Determining the turbulent lifetime then 
amounts to computing the mean first passage time at the 



absorbing boundary e = £ m i n - 

Such a calculation is however difficult for an arbitrary 
functional form of the entropy s(e, R) and we have to 
restrict the choice of s{e,R) to the linear form 

s(e, R) = MR) (e - e mi „(i?)) + s (R) (12) 

over the interval £ m i n (i?) < £ < s max (R). In this case, 
the mean first passage time can be computed from the 
solution of the associated Fokker-Planck equation [75], 
and one finds 



V 



(Ae) 2 f(v(f3g-f3)Ae) (13) 



with Ae — £ max — £min and fi g — j3{R g ), and where the 
function fix) is given by 



1 



1 



(14) 



For large V, the argument of the function / in Eq. (131 is 
large as soon as P g ^ P, that is R ^ R g . The value of r is 
then given, to a good approximation, by the asymptotic 
behavior of fix) when x — > ±oo, which reads 



fix) 
fix) 



1 

R 

e x 



x — > — oo , 
+00 . 



Hence, r is given for /3 g < (3 by 

A£ 



D {P-P a 



(15) 
(16) 

(17) 



which turns out to be independent of the volume V , as 
intuitively expected in the large V limit. In terms of 
Reynolds number, one thus has a power-law divergence 
for R close to R g [R < R g ), 



Rn — R 



(18) 



However, for any finite volume V this divergence is cut 
off when R approaches R g , as soon as R g — R < aV^ 1 
with some constant a, and a crossover is observed to the 



exponential form obtained from Eq. (16) 



DoViJlg ~ PY 



(19) 



Contrary to Eq. ( 17 1, expression ( 19 1 involves the volume 
V . For V — > oo, r becomes infinite, and a true power-law 
divergence is observed for R < R g . For very large but 
finite V, the divergence can be observed in practice only 
on a narrow range of Reynolds number, before r becomes 
exceedingly large. On this narrow range, (/3 g — P)Ae 
behaves linearly with R. In contrast, if V is not too 
large, the range of R over which the divergence is ob- 
served broadens, and corrections to the linear behavior 
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FIG. 4: (Color online) Sketch of the turbulent lifetime as a 
function of /3 g /j3, the effective Reynolds number. Main panel: 
the continuous curve is the lifetime for a given volume of the 
system as given by Eq. | |13[ ); dashed curves are the asymptotic 
functional forms which govern the behavior of r on each side 
of P g -they have been shifted for clarity. Inset: Turbulent 
lifetimes for increasing system size: the larger the system, the 
steeper is the increase of lifetime. The singular behavior is 
observed in the limit of infinite system size only. 



of (/3o — |5)Ae with R can become observable, possibly 
leading to a super-exponential behavior of r as a func- 
tion of R. Though sub-exponential behavior cannot be 
discarded, one expects at least Ae to increase with R, 
which goes in favor of the super-exponential case. 



IV. DISCUSSION 

What did we learn from the physics of glasses? 

The initial motivation of the analogy proposed in this 
paper was two-fold. First the intense debate that ani- 
mated the transition to turbulence community regarding 
the possible divergence of the turbulent lifetimes at a 
finite Reynolds number was reminiscent of a similar situ- 
ation encountered in the physics of glasses a few decades 
earlier. Second, the idea that the dynamics is controlled 
by unstable solutions away from the laminar state shared 
some similarity with the role played by the large number 
of saddles at the onset of the glass transition. The goal of 
the analogy presented here was to make these intuitions 
more precise. 

We have shown that indeed, even with very good 
data, one cannot discriminate a singular dependence from 
a regular but very fast increase of the turbulent life- 
times, especially if one includes the possibility of a Vogel- 



Fulcher-Tammann like singularity. We have also seen 
that finite size effects may lead to a crossover, which can- 
not be resolved experimentally or numerically because of 
the extremely large timescales at play. It is remarkable 
from that point of view that the model we introduced ex- 
hibits a crossover between the two main functional forms 
debated in the literature, namely the power-law and the 
super-exponential ones. 

The model presented here was designed to be as simple 
as possible, and for the purpose of illustrating the anal- 
ogy. As such, it does not pretend to be realistic in any 
way. Valuable improvements could be obtained by mea- 
suring in direct numerical simulations the dissipation rate 
as a function of the energy density, as well as characteriz- 
ing the statistical properties of the turbulent energy fluc- 
tuations in the intermediate range of Reynolds numbers. 
In addition, counting the number of unstable solutions 
as a function of their energy density, that is accessing 
s(e, R), would be a major step towards the characteri- 
zation of the transition to turbulence. This is obviously 
a difficult task, but still far less ambitious than char- 
acterizing the stability properties of these solutions and 
describing the complex interplay of their stable and un- 
stable manifold. This simplification is in essence the one 
gained when switching from a dynamical systems point 
of view to a statistical physics one. 

In the above section, we have considered V as the vol- 
ume of the system. However, in the spirit of real space ap- 
proaches the relevant volume to be considered may rather 
be the volume of coherent regions of the flow, namely re- 
gions over which correlations extend. In a very large 
aspect ratio experiment, it is plausible (though not ob- 
vious) that far away regions in the system experience no 
interactions. As a result, the volume V would acquire a 
more intrinsic nature: it would then be self-determined 
by the flow dynamics and not by the arbitrary size of the 
experiments. 

Such a coherence volume cannot be accessed in the 
framework of models similar to the Random Energy 
Model, which is mean-field in nature. However, if the 
analogy with the physics of glasses proves to be fruitful, 
it would be of interest to consider its most recent devel- 
opments (including in particular the Random First Or- 
der Transition scenario [76 J, which precisely address the 
real-space description issue [23] • Pomeau [15] suggested 
more than twenty-five years ago that the growth and 
death of the laminar and turbulent regions could obey 
a first-order nucleation-like dynamics (albeit of a pecu- 
liar type, given the fluctuating active property of the tur- 
bulent state and the adsorbing character of the laminar 
state). Let us conclude with the somewhat naive sugges- 
tion that taking inspiration from the Random First Or- 
der Transition theory of glasses might be a way to extend 
the standard laminar-turbulent coexistence scenario to a 
situation where a large number of turbulent states (as- 
sociated to local unstable solutions of the Navier-Stokes 
equation) coexist with the laminar state. 
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